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We analyze relationships between quantum computation and a family of generalizations 
of the Jones polynomial. Extending recent work by Aharonov et al., we give efficient 
quantum circuits for implementing the unitary Jones- Wenzl representations of the braid 
group. We use these to provide new quantum algorithms for approximately evaluating 
a family of specializations of the HOMFLYPT two- variable polynomial of trace closures 
of braids. We also give algorithms for approximating the Jones polynomial of a general 
class of closures of braids at roots of unity. Next we provide a self-contained proof of a 
result of Freedman et al. that any quantum computation can be replaced by an additive 
approximation of the Jones polynomial, evaluated at almost any primitive root of unity. 
Our proof encodes two-qubit unitaries into the rectangular representation of the eight- 
strand braid group. We then give QCMA-complete and PSPACE-complete problems 
which are based on braids. We conclude with direct proofs that evaluating the Jones 
polynomial of the plat closure at most primitive roots of unity is a #P-hard problem, 
while learning its most significant bit is PP-hard, circumventing the usual route through 
the Tutte polynomial and graph coloring. 
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1 Introduction 

There is evidence that a computer which could manipulate quantum mechanical degrees of 
freedom would be more powerful than a classical computer. Since the discovery by Shor that 
a quantum computer could efficiently factor composite integers, a task which is believed to 
be hard on conventional computers, much effort has been expended toward understanding 
the capabilities and limitations of quantum computers. In this paper, we investigate ways 
in which the Jones polynomial invariant of knots and links, together with its generalizations, 
contribute to this understanding. Wittcn discovered |40] that the Jones polynomial could 
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be understood via tools from topological quantum field theories. Freedman, Kitaev, Larsen 
and Wang [11] established its connection to computer science in the context of topological 
quantum computation. In a topological quantum computer, the trajectories of particles which 
are restricted to a plane are braided in order to manipulate the internal state of the computer. 
They showed that such computers can be simulated on conventional quantum computers, and 
also that when the underlying physics of a computer is described by a suitable topological 
quantum field theory, that it is possible to achieve universal quantum computation. 

In [lOj . it was shown that topological quantum computers could be efficiently simulated 
by computers based on the standard quantum circuit model, implying the existence of an 
efficient algorithm for approximating the Jones polynomial of certain links obtained from 
braids. More recently, Aharonov, Jones and Landau [2] gave explicit quantum algorithms for 
approximating the Jones polynomial of either the trace or plat closure of a braid on a quantum 
computer. Our first contribution in this paper is a family of quantum circuits implementing 
the so-called Jones- Wenzl unitary representations of the braid group, generalizing the circuits 
for the Jones representation from [2]. We then show how to use our circuits to approximate 
a family of evaluations of the HOMFLYPT two-variable polynomial. We also sketch two 
different quantum algorithms for approximating the Jones polynomial of links obtained via 
a class of closures generalizing the trace and plat closures. In the sense of j7j, our quantum 
algorithms obtain additive approximations of Jones and HOMFLYPT polynomial evaluations 
in polynomial-time. 

In [121 [13], ^ converse to the above simulation results was proved. More specifically, it was 
shown there that quantum circuits can be simulated by braids, in such a way that the output 
probabilities of the quantum circuit are functions of the Jones polynomial of the plat closure 
of some braid of comparable length to that of the original circuit. Our third contribution is a 
simpler proof of their main result. In fact, this work grew from our attempts to understand 
the connections between these works and [2]. 

Finally, we give applications of the Jones polynomial to quantum complexity theory. We 
begin by restating a result that was proved in [J , showing that a machine which obtains an 
additive approximation of the Jones polynomial of the plat closure of a braid is equally as 
powerful as a quantum computer. We then introduce two new problems, Increase Jones 
Plat and Approximate Concatenated Jones Plat. The first asks if a given braid can 
be conjugated by another braid from a given class such that the Jones polynomial of its plat 
closure is nearly maximal. We prove that this problem is complete for the complexity class 
QCMA, a certain quantum analog of NP. The latter problem asks if a given braid, after being 
concatenated with itself exponentially many times, has a large plat closure. We show that 
this problem is PSPACE-complete. Finally, we give self-contained proofs that learning n of 
the most significant bits of the Jones polynomial of the plat closure of an n-strand braid is 
a #P-hard problem, while learning its most significant bit is PP-hard. This constitutes a 
simpler, quantum-based proof of the known result that computing the Jones polynomial is 
#P-hard. The original proof of this (see e.g. [39]) relates the Jones polynomial to the Tutte 
polynomial of a signed graph to the chromatic number of that graph to the #P-complete 
problem #SAT. Our proof, however, shows that learning a linear number of the highest order 
bits of the Jones polynomial of the plat closure of a braid is ^^P-hard. 

The paper is organized as follows. In the next section, we review elements of the theory 
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of links, braids and the Jones polynomial. There, we review some different ways of turning 
a braid into a link, while introducing a new type of elosure which we call a generalized 
closure. We also give precise statements of the approximations achieved by the algorithms. 
In Section [3l we describe the unitary representations of the braid group used in this paper, 
highlighting their origins from Hecke algebras. In Section IH we show how to efficiently 
implement these representations on the state space of a quantum computer. Section[5]contains 
derivations of representation-theoretic formulae for the HOMFLYPT and Jones polynomials. 
Section [6] shows how to use the quantum circuits of the previous section to construct quantum 
algorithms for approximating certain evaluations of the HOMFLYPT and Jones polynomials 
of closed braids. In Section [71 we provide a self-contained proof of the converse of the results 
in the previous section - that local quantum circuits can be simulated by braids. Finally, 
in Section [SI we review notions of classical and quantum complexity theory, after which we 
present our complexity-theoretic results. We conclude with a discussion in Section [5) 

2 Background on link invariants and additive approximations 

Here we give background material on the links invariants used in this paper. We also formalize 
the notion of approximation for our quantum approximation algorithms. We refer the reader 
to. for example, |39| for further details on link invariants. A knot is a closed, nonintersecting 
curve in M.^. More generally, a link is an embedding of some finite number of nonintersecting 
closed curves in R'^. Links are identified up to isotopy, which means that two links which 
arc related by some bijection / : M'^ IR'^ of the ambient space with itself for which / and 
are continuous are considered to be equivalent. An oriented link is a link in which every 
component is assigned an orientation. If L is an oriented link, we will denote the unoriented 
link resulting from forgetting the orientation as L. 

A central problem in knot theory is to determine, given descriptions of two links, whether 
or not they are equivalent. In order to solve such a problem, one must decide on how the 
links are to be described. One such way is in terms of a link projection. Informally, this 
is a two-dimensional diagram which uniquely specifies the link up to isotopy. The reader is 
referred to textbooks such as [25l [28l [29l [39] for examples and further exposition. The key 
property of any link projection is that at each crossing, it keeps track of which string goes 
above or below the other. Usually, this information is conveyed by leaving a gap in the string 
of the undercrossing. A classic result in knot theory states that two unoriented links are 
equivalent if and only if the link projection of one can be transformed into that of the other 
by a finite sequence of Reidemeister moves (see e.g. any of the texts listed above). 

The HOMFLYPT and Jones polynomials 

Another way of (partially) distinguishing links involves assigning an invariant to each link 
so that links which are equivalent up to isotopy have the same invariant. The HOMFLYPT 
polynomial [16i 132] Hj*{t,x) G Z[t*^,a;*^] is one such invariant. It is inductively defined by 
the skein relation 

tHL+{t,x) -t-^HL_it.x) ^ xHL„{t,x) (2.1) 

(see Figure [l?!]) together with the fact that the HOMFLYPT of the trivial knot (that is, a 
single unknotted circle) is unity. In this paper, we are concerned with various one-variable 
specializations of this invariant. For each integer k > 2 there is a corresponding one-variable 
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Fig. 2.1. L_|_, L_, Lq denote three oriented links that differ in a smaU region, as symbohzed by 
the above diagrams, but are otherwise the same. We say that the crossing corresponding to L+ is 
a positive crossing, and the one corresponding to L_ is a negative crossing. 



HOMFLYPT polynomial in Z[q±i/2]: 

i/f (g)=i/^(,^-/2,^i/2_g-i/2) (2.2) 

This is also known as the s[fe invariant. Restricting to fc = 2 yields the famous Jones polyno- 
mial: 

If L consists of m unlinked copies of the trivial knot, setting t = q*^/^ and x = q^^^ — q^^^^ 
in (|2.ip implies that 

rt.^ /„fe/2 „ '"-I 

The exponentiated quantity on the right is known as a quantum integer and is easily shown 
to satisfy 

qi/2_l-i/2 =1^+ 1^ (2.3) 

a q ^ 1. By continuity, we consider this expression to hold as q ^ 1, recovering the ordinary 
integer k. In this paper, we will only be interested in evaluations of these invariants at 
primitive roots of unity q = e^'^'/^, where i is a, nonnegative integer. To this end, we define 
for each integer k and each nonnegative integer £, the constani|° 



_ gfc/2 _ q fc/2 



sin{nk/i) 
g^e^.^/i sin(7r/^) ■ 



A value which will occur frequently, especially when analyzing the Jones polynomial, is 

[2], = 9^/2 + 9-1/2 =2cos(7r/^). (2.4) 

The Jones polynomial will be the central object of our study, although we will also provide 
algorithms for approximating the one-variable HOMFLYPT polynomials. Where convenient, 
we sometimes write J( L , 9) = J-^{q). 



'^Our definition here departs from the usual convention 1381 1121 which defines \k\q = 1/2 — -1/2 ^^'^ arbitrary 
q. This expression evaluates to ours when g is a primitive £'th root of unity, which is the case of interest to 
us. We warn that other authors (e.g. 1191 ) use a different convention for the quantum integer and q. 
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The braid group 

The braid group i?„ on n strands is generated by eounterclockwise twists {a'i,tT2, 
whieh satisfy the relations 



(JiCTj = fT.jCr,; 



(2.5) 
(2.6) 



The reader should picture n hanging pieces of string, and interpret each generator fj i as a 
counterclockwise exchange of the strands in positions i and i+l, with inverses of the generators 
corresponding to clockwise twists. For instance, the braids ai,(72^ ^ ^3 given by 



and 



r 



The group product then corresponds to vertical concatenation of braids (up to isotopy equiv- 
alence). For instance, 



(Tl = 



The relation (|2.5p is known as the Yang-Baxter equation, which is diagrammatically expressed 
as 



\ 

For each n, i?„ is isomorphic to the subgroup of Sn+i generated by {cri, . . . , (Jn-i}, consisting 
of braids in which the n + I'st strand does not participate. We denote the associated inclusion 
map as t: i?„ Sn+i, which acts by adding an extra strand to the right of any given braid. 
Given a braid h € Bn, there are many possible ways of closing the braid b to obtain a link. 
In the following we will recall three possible ways, the trace, the plat and a combination of 
these we call a generalized closure. 

Specifying links with braids 

One way to turn a braid into a link takes, for each 1 < i < n, the i'th strand at the top of 
the braid and glues it to the i'th strand at the bottom. The resulting link h is called the trace 
closure of h. For instance, the trace closure of the braid b = (72(7^ ^(72 is isotopic to the union 
of the trivial knot and the Hopf link: 
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Every link can be represented as the trace closure of some braid. Conventionally, one considers 
the trace closure of a braid to yield an oriented link, since an orientation can be unambiguously 
defined to travel downward along the braid. With this convention every generator ai gives rise 
to a positive crossing, while its inverse cr"^ to a negative crossing in the sense of Figure EH] 
A theorem of Markov demonstrates that, given braids b and b', their trace closures b and 
b' are isotopic if and only if b and b' are related by a finite sequence of Markov moves. We 
review these moves in Section [5] as they are relevant for obtaining a representation-theoretic 
formulae for the one- variable HOMFLYPT and Jones polynomials of trace closures of braids. 

Another way of turning a braid on an even number 2p of strands into a link is known as 
the plat closure. Such a procedure starts with a braid b G i?2p and connects adjacent pairs of 
strands at the top and at the bottom. For example, the plat closure of = a2a'^^a'2 & B4 is 
isotopic to the trefoil knot, as 




Given a braid b € Bn , we define a generalized closure of b to be specified by the following 
data: two braids x,y € i3„ and nonnegative integers p and r which satisfy 2p + r = n. The 
generalized closure Xx'yi^) of ^ obtained by first forming the braid xby, performing the plat 
closure on the 2p leftmost strands at the top and bottom, then connecting the remaining r 
strands at the top with the corresponding r strands at the bottom as in the trace closure. 

Unlike with the trace closure, there is no unambiguous way to assign an orientation to the 
plat or generalized closure of a braid (unless n — t). While this poses a potential problem for 
defining the Jones and HOMFLYPT polynomials of such closures of braids, it is known [55] 
(and easily proved using the Kauffmann bracket [21] formula) that if L and L' are two oriented 
links that are isotopic as unoriented links, their Jones polynomials satisfy J-j^i<l) = Q"^Jj^i<l) 
for some easily computable m G Z. This implies the absolute value of the Jones polynomial at 
roots of unity is independent of the orientation of a link, and is thus an invariant of unoriented 
links. In Section [5] we provide a representation-theoretic formula for the absolute value of 
the Jones polynomial of any generalized closure of a braid at primitive roots of unity, while in 
Section [6] we present a quantum algorithm which approximates these invariants of unoriented 
links. 

Additive approximations 

The notion of approximation which we will require was recently formalized in [7] ; we review 
that material here. Let A" be a set of problem instances, and suppose we are given a non- 
negative function f : X which is potentially difficult to evaluate exactly. The main idea 
of [3 is to approximate the function / with respect to some positive normalization function 
g: X ^ W. Departing slightly from the definition given in [7| , an additive approximation 
for the normalized function f /g associates a random variable Z{x) to any problem instance 
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X ^ X and 6 > satisfying 



Pr 



-Z{x) 



<5\> 3/4. 



In addition, it is required that such an approximation be achieved in time which is polynomial 
in the size of the problem instance and in 1/5. 

In Section [5] we provide tight upper bounds on the evaluations of the invariants we wish 
to approximate. We will see in (|5.24p that li h e i?„, the one- variable HOMFLYPT of the 
trace closure obeys |iji''^(e2"/^)| < [fc]""^ Similarly, we show in ([OS]) that ii n ^ 2p + t 

and x,y G the Jones of the generalized closure satisfies \J{£,x'Jy{h),e'^'"'^/^)\ < [2]^"*"*"^. 
These bounds are achieved by closing the identity braid for the one-variable HOMFLYPT 
of the trace closure, or otherwise the braid x~^y~^ for the Jones of the generalized closure. 
We use these normalizations for our additive approximations, so that the absolute values of 
the resulting normalized polynomials always lie between and 1, regardless of the size of the 
braid, number of strands, or particular primitive root of unity. We therefore consider the 
following two problems: 

Problem 2.1 (Approximate HOMFLYPT Trace Closure) Given is a braid b £ Bn of 
length m, S > and positive integers fc, £ satisfying 2 < k < (., where k is a constant. The 
task is to sample from a random variable Z £ C, \Z\ < 1, which is an additive approximation 
of the one-variable HOMFLYPT polynomial of the trace closure, evaluated at e^'^'/^, in the 
sense that 

1 



Pr 



-Hi'\e 

h 



2iri/ 



- z 



<S}> 3/4. 



Problem 2.2 (Approximate Jones Closure) Given is a braid b E _B„ of length m, two 
braids x,y G _B„ of length 0(poly(TO)), positive integers £,p and r which satisfy 2p + r = 
n, and S > 0. The task is to sample from a random variable < Z < 1 which is an 
additive approximation of the absolute value of the Jones polynomial of the generalized closure, 
evaluated at ^'^'<-l^ in the sense that 



Pr. 



1 



[2] 



p+r — 1 



- Z 



<&)> 3/4. 



In Section [SI we prove: 

Theorem 2.1 Each of the above two problems can be solved in 0(poly(m, 1/5)) time on a 
quantum computer. 



3 Unitary Jones-Wenzl representations of the braid group 

Here we summarize the representation theory of the braid group that we will need, following 
[38j . The relevant representations are inherited from representations of a certain family of 
Hecke algebras Hn{q) which generalize the group algebra C5„ of the symmetric group 5„. In 
Section [51 we will show how these representations give rise to useful formulae for evaluations 
of the Jones and one-variable HOMFLYPT polynomials at roots of unity. 
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Fig. 3.1. Young diagram for the partition A = [5,3,2]. 



The Hecke algebra 

It is well-known that the group algebra C5'„ of the symmetric group S'„ on {1, 2, . . . , n} objects 
has a presentation in terms of generators {1, si, S2, ■ ■ ■ , s„_i}, where s,; is the involution which 
swaps i and i + 1. These satisfy the relations 

= 1 (3.1) 

SiSjSi = SjS^Sj, \i-j\ = l (3.2) 
s,Sj = s-jSi, \i - j\ > 1. (3.3) 

For every q G C^, the Hecke algebra Hn{q) of type A„_i is defined to be the algebra over C 
generated by 1 and {gi, . . . ,gn-i} satisfying the relations 

gf = g^{q-l) + q (3.4) 
9t9i9t = 939t9j = l (3-5) 

9i93 = 9]9t j\ > 1- (3-6) 

Note that this is a deformation of the group algebra CSn of the symmetric group Sn, which 
is obtained when q = I. While it is known that H„{q) ~ CSn whenever q is not a root of 
unity, we will rather be interested in the cases when g is a primitive root of unity, where the 
representation theory is slightly more subtle than that of CSn- One may represent the braid 
group Bn inside the Hecke algebra Hn{q) by simply mapping each generator (t.; of i?„ to the 
generator gi of Hn{q). Indeed, the representations of _B„ we will use in this paper are induced 
by representations of i?„(g) on suitable finitc-dimensional Hilbcrt spaces. As with the braid 
groups, for each n, Hn{q) is isomorphic to the subalgcbra of iJ„+i(g) generated by 1 and 
{gi, . . . , g„_i}. We use the same notation for the inclusion maps as with the braid groups, 
writing l: Hn{q) — *■ Hn+i{q)- One important distinction is that, writing 1„ for the identity 
in Hn{q), we have i(l„) = In+i- We will therefore often omit the subscript for the identity 
element of each Hn{q) whenever there is no cause for confusion. 

Young diagrams and tableaux 

As with Sm the unitary irreps of _B„ are labeled by partitions A = [Ai, A2, . . . , A„] satisfying 
Ai > A2 > ■ ■ ■ > A„ > and Ai = n. Throughout, we consider all partitions obtained by 
adding or deleting trailing zeros to be equivalent. We will identify any such partition with its 
Young diagram as pictured in Figure [5TT] for the partition A = [5, 3, 2]. Note that the diagram 
has one row for each part of the partition, and that the i'th row contains A^ boxes. We allow 
for the empty partition and diagram, which we denote 0. Let A„ be the collection of n-box 
Young diagrams. Given two diagrams A € A„ and € A™, their sum A-I-/X G J^n+m is always 
well-defined, having parts equal to (A + /i).; = A.; + ^i. When it is well-defined as a partition, 
we will similarly define the difference A — /i G A„_,„ to have parts (A — ~ Xi — fii. 
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□ 




Fig. 3.2. The Young graph Ao,4. 



Define the sets of diagrams Am,n = 1+1."=™ and abbreviate A = Ao,oo- We may consider 
elements of A to be the nodes of an infinite directed graph, which we will call the Young 
graph, which contains an edge from a diagram A' to a diagram A whenever a single box can 
be added to A' to obtain A. We remark that we depart slightly from common terminology, 
where this graph, with each arrow reversed, is called a Bratteli diagram. With a slight abuse 
of notation, we write A for both the set of diagrams and for the corresponding Young graph, 
making a similar identification between subsets of A such as Ao,„ and the associated subgraph. 
In Figure [3^ we give a picture of Ao,4. 

By a numbering of an n-box diagram, we mean any assignment of the numbers {1, 2, . . . , n} 
to the boxes of the diagram, where each number appears only once. For a given diagram 
A G A„, we write T\ for the set of standard tableaux with shape A, corresponding to numberings 
which are strictly increasing along each row and column. We often call a standard tableau 
just a tableau. We denote the set of all n-box tableaux as r„ = l+J;^^^ T\ and the set of 
all tableaux as T = r„. We identify members of T„ with length-n paths starting at 

in Ao,n, pictured in Figure [3721 with elements of T\ correspond to paths which end at the 
diagram A. 

In this paper, we shall rather require certain restricted classes of Young diagrams and 
tableaux. To each pair (/c, t) of integers satisfying £ > fc > corresponds a class of irreducible 
unitary representations of These representations are labelled by the 7i-box [k, £)- Young 
diagrams, defined as 

Al'^^^) = {A e A„ : Xk+i =0,Xi-Xk<t- k}. 

Note that these are the diagrams with at most k rows for which the difference between the 
numbers of boxes in the first and fc'th rows is at most £ — k. We will often just say that such a 
A is an (k, ^)-diagram. For a given diagram A, we refer to Ai — Xk as the level of A, so that i—k 
is the maximum level of all diagrams in A^*^'^^. As before, we make the similar abbreviations 
Am.'n and A'*^'^^. We also identify any of these sets of diagrams with the appropriate subgraph 
of A. We refer to such a graph as a (fc, £)- Young graph. In Figure 13. 3[ wc show the graph 
Ag^l^-*. For a given A G A^i'^\ define the {k, £) -tableaux T^'^'^ C T\ to be those standard 
tableaux of shape A for which successively deleting the largest numbered boxes yields, at each 
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Fig. 3.3. The Young graph Comparison with Figure [3^21 reveals that this is the subgraph 

of Aq 4 obtained by removing vertices whose diagrams either have more than two rows, or level 
greater than 3 = 5 — 2. 

step, another (fc, ^)-diagram. These tableaux can be identified with paths in 1^^^!^ going fi'om 
to A. We also define T^t'^^ and T^*^'^) in a similar manner as above. In the Section [5l we will 
see that the (2, ^)-diagrams and tableaux are what is needed to give representation-theoretic 
expressions for the Jones polynomials J{e^^^^^) of closures of braids, while the (A:, £)-diagrams 
are relevant for the one-variable HOMFLYPT polynomials H^'^^ [e^^^^^). 

Unitary Jones- Wenzl representations of Bn 

The representations of _B„ we will introduce arc parameterized by integers £ > A: > 0, where 
£ specifics a primitive .^'th root of unity which we will denote q — g^'^*/^. The irreducible 
components of the representations we will describe reduce to the usual irreps of Sn when 
k = n and £ oo, so that q 1. For the rest of this section, we fix £ > fc > and 
set q = Q^TTi/i^ describing a corresponding representation of Hn{q) which induces a unitary 
representation of Bn- This representation is more easily described by a change of variables 
which leads to an equivalent presentation of i/„ (q) . Rewriting the quadratic relation (|3.4p as 
(<7i + ^9) = 0, we may define idcmpotents e; = (9 — <7i)/(l + 9) (meaning that ef = Ci) 
for which 

gi = q{l - Ci) ~ ei ^ q - (1 + q)ei. (3.7) 

Using these idempotents, the relations p.4H3.6p may be expressed as 

e,? = e, (3.8) 
616^-61 TCi ejCiBj - TCj, \i- j\ = I (3.9) 

BiCj = Cj-ej, |z-j|>l, (3.10) 

where wc follow the usual convention in setting r = [2]^^. For each integer d, define 

«M = (3.11) 

noting that ai{d) + ai{—d) = 1 = (1). Given any n-box tableau t and an integer 1 < i < n, 
we respectively write Ci{t) and ri{t) for the column and row which contain the number i, 
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defining 

d^it) = c,{t) - c,+i{t) - {n{t) - n+i{t)). (3.12) 

Note that di{t) gives the total number of leftward and downward steps it takes to "walk" 
from i to i + 1 on t, where moving right or up counts negatively. Given a generator Si of Sn 
and a (A;, £)-tableaux T^'^\ let t' be the numbering obtained by swapping the numbers i and 
i + 1 vat £ t!^'^\ If t' is also a (fc, £)-tableau, define Siit) = t' and Si{t') = t; otherwise set 
Si{t) = t. In order that Siit) ^ t, it is necessary and sufficient that i and i + 1 not be in the 
same row or column of t. In such a case, we may think of it as a deformation of the path in 
Ag ^ corresponding to t which changes only the i-box diagram passed in the original path. 

(k i) . 

Note that each generator Si induces a partition of T„ into sets of numberings (or paths) of 
size at most two. 

(k £) / \ (k £) 

We may now define a representation 7r„ of i7„ [q) on a vector space Vn ' which has an 
orthonormal basis labeled by n-box (fc,£) tableaux, i.e. 

y^^)=span{|t):ier^^)}. 

For each i, tt^'^'^-' (e^) is a projection which is block diagonal in the T^'^'^' basis with blocks of 
size either one or two. The size two blocks correspond to distinct pairs of tableaux t and Si{t) 
which are contained in Tn'''^\ Letting Vi^t = span{|t), \si(t))}, the restriction of 7ri'^'^'(ei) to 

Vi^t is the rank one projection onto ai{di{t)^\t) + ac(^—di{t)) \si{t)), written 




a({d,it)) ^ae{d,{t))ae{-d,{t)y 
ae{dt{t))ae(~di{t)) ai[-di{t)) 



(3.13) 



For the remaining i e r„ ' , set 7r„ ' ' {ei)\t) = \t) if z and i + 1 are in the same column of 
t; otherwise set 'Kn'^\ei)\t) = 0. This defines the action of Tr^'^\ei) on all of Vn'''^\ Using 
p.7p . this representation maps each gi to the unitary matrix 

4'=.^) (5.) = 9V<..) - (1 + q)ni'''He^)- (3.14) 
Abbreviating d — di{t), we may use (|2.3p to express the restriction to each Vi^t as 

-it^'Hg.)^ . = -i'^' ( ^^^^^ (3-15) 

If i and i + 1 are in the same column of this gives TT'^'^\gi)\t) = —\t), while if they are in 
the same row we have 7rl'^'^^((7i)|f) = q\t). Because permutations cannot change the shapes 
of tableaux, it follows that this representation is reducible. The irreducible components are 
labeled by the (A:, ^)-diagramS; so that 



where each n-^ ' acts nondegenerately only on the subspace ' ~ {\t) : t £ ' }, and 



AeA 



{k,i) 

A 
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4 Efficient quantum circuits for unitary Jones- Wenzl representations 

We will now show, for any integers I < k < £, how to implement the corresponding Jones- 
Wenzl representations tt^'^'^'' defined in (j3.14p on a quantum computer. The basic idea is to 
embed these representations into a tensor product space (c.f. Section 5.2 of [33]). Our circuits 
resemble those of [2], where the k = 2 case is treated, although we require the combinatorial 
tools of Young diagrams and tableaux introduced in the previous section for higher values of 
k. These circuits allow us to give a new quantum algorithm in Section [6] for approximating 
the family of one-variable HOMFLYPT polynomials TJ^*^) of the trace closures of braids. We 
begin with some background on quantum circuits. 

The state space of a quantum computer has the structure of a Hilbert space Ti.. Denoting 
by UiTi) the group of unitary matrices acting on 7i, a quantum gate is any U € U(Ti.). By 
the standard quantum circuit model, we shall mean quantum computers whose state space 
decomposes into a finite number of localized two-dimensional subsystems, or qubits, so that 
H = (C^)*^". We fix a preferred orthonormal basis for the Hilbert space of each qubit, which 
we call the computational basis: 




The state space of Ti. is thus spanned by tensor products of these basis vectors, which we 
write as {|a;") : x" G {0, 1}"}, where 

Ix") EE |a;i)|a;2) • • • |a;„) ee |a;i) (g> \x2) «)•••«) |a;„). 

We will abbreviate the set of all gates on n qubits as ?7(2") = J7((C^)®"). We assume that 
the qubits are arranged on a line so that each qubit has at most two neighbors. In this case, 
it is well-known that an arbitrary unitary U e J7(2") can be performed on the qubits by a 
sequence of (possibly exponentially many in n) local unitaries from U{4) which act only on 
adjacent pairs of qubits. By a length-m quantum circuit, we will mean a sequence of m local 
two-qubit gates. When appropriate, we identify a quantum circuit with its overall unitary 
transformation U ~ UmUm-i ■ ■ - Ui, where each Ui is a two-qubit unitary. Below, we will 
work with quantum circuits whose local systems (registers) may have dimension k > 2. Such 
systems can be implemented using O(logfc) qubits, on which any unitary can be performed 
using only 0(k) two-qubit gates. 

We now proceed by describing an encoding of the basis states of the reducible represen- 
tation FC^''^) into 0(n log(fc)) qubits. Then we show how, using 0(nlog(n)) ancillary qubits, 
one can implement the image Tr^'^\ai) of each braid group generator ai with a quantum 
circuit of length 0(poly(n, fc)) . Given a tableau t € T^'^'^'' and an integer I < i < n, recall 
that we write ri{t) and Ci(t) for the row and column which respectively contain the number i. 
Notice that these numbers satisfy 1 < ri{t) < k and 1 < Ci{t) < n for each i. We abbreviate 
r{t) ~ ri{t)r2{t) ■ ■ ■ rn{t) £ {1, . . . , fc}" and define c{t) G {1, . . . , n}" similarly. Our first ob- 
servation is that r{t) uniquely specifies t. We may thus assign to each t a computational basis 
state of n systems with Hilbert spaces 7?,i, 7?,2, . . . , TZn, where each Hilbert space TZi ~ C'^ is 
spanned by computational basis states {|1)^% . . . , |fc)^*}. Together, the systems have a com- 
bined Hilbert space which we write 7^" (C'')'^". This yields an embedding \t) ^ \r{t))'^" 
of Vn'^'^^ into TZ" whose image we call the computational subspace. 
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We now introduce, for each l<i<77, — l,a register with Hilbert space Vi which will 
store the value di{t) defined in (|3.12p for each tableau. Because < n, we let Vi be 

spanned by basis states {| — n)^% \ —n+ 1)^', . . . , We may now define an isometry 

W : 7?," 7?."I?"~^ which coherently computes these values for a given encoding of a tableau 
as 

Wlrit))-^" ^\r{t))'^"\d{t)f"-\ (4.16) 

Below, we will show that this isometry can be implemented efficiently. Now we follow (|3.15p to 
define the unitaries which perform the braiding, while taking care to appropriately update the 
registers. For each 2 < i < n—2, define the unitary Ul which acts on TZiR-i+iDi^iDiDi+i 

as 



where aa = — ^ ^^j^ and (3d = ^Q^^"^ ~ j^- Since di can only equal ±1, we also set 



u[\d,r^ 



ifdi = l 



Additionally, we define J7^_i to act on 7?.„_i7?.„P„_22-'n-i via 

Each of these unitaries can be implemented using 0(poly(n, fc)) two-qubit gates. Define 
the unitaries Ui = W~^U^W for 1 < i < ti. By construction, the map i-^ Ui defines a 
(degenerate) unitary representation pil^'^^ of i?„ which is isomorphic to Trll''^\ Suppose now 
that we are given a braid b = ■ ■ ■ o-^^o-^^^ G Bn- The corresponding unitary 

X\ TjXrrt TTX2TTX-\ 

therefore evolves the computational subspace of 7?." according to the following commutative 
diagram 

\ t 

-R" — ^7^". 

Now we show that the isometry W in (|4.16|) can be efficiently implemented. For this, we 
introduce n systems with joint Hilbert space C" , where each Ci is spanned by computational 
basis states {|0)'^', . . . , \n)'^'}. The purpose of these systems is to store the column numbers 
of each i. In addition, we introduce k registers £1, . . . , Ck with bases {|0), . . . , \n)} which will 
store the shape of the tableau. Given a qubit encoding |r(t))^ of a tableau t, we define the 
isometry W : Tl" WC'C^ as 

W^'|r(t))^" = |r(t))^'>(t))'^"|A(t))^\ 
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where X{t) is the vector of row lengths of t. This isomctry can be efficiently implemented 
by first preparing the additional systems in the state |0---0)'' |0---0)''' , after which the 
composition of local gates W' — VF,' • • • is applied, where each W[ acts on TZiCiL^ as 

where A' is A with an additional box in row r. For each \ < i < we then define local 
unitaries Wi which act as TZiTZi+iCid+iVi as 

W.lr)"^^ |c)C. |c')C'+i \0f' = |r)^' \r')^^+^ \cf' \c'f'^' \c-c'-{r- r'))^' . 

We may thus write W = W'^^Wn^i ■ ■ ■W2W1W' , where we uncomputc the C" and 
registers. Because each Wi and can be implemented with 0(poly(n, fc)), we find that W 
can be implemented with 0(poly(n, fc)) local gates. 

5 Representation-theoretic formulae for the Jones and HOMFLYPT 

We will now illustrate how to use the unitary representations of Section [3] to write expressions 
for the evaluations, at primitive roots of unity, of the Jones and one-variable HOMFLYPT 
polynomials of the trace closure of a braid. We then review the Temperley-Lieb algebra, which 
we use to derive a formula for the Jones polynomial of our generalized closure of a braid. 

Jones and HOMFLYPT polynomials of trace closures from the Markov trace 

Here, we review the theory of Markov traces on Hecke algebras and show how they yield 
formulae for the one-variable HOMFLYPT polynomials of the trace closure of a braid. Our 
presentation is primarily based on Section 6 of [19] and Section 3 of [38] . We urge the 
reader who seeks a deeper understanding to consult these references. We also recommend [4]. 
Throughout this section, we fix a primitive root of unity q = e^"^^!^ for some integer ^ > 3. 

The theory of Markov traces is best understood in light of the so-called Markov moves on 
braids. As mentioned earlier, the trace closures of two braids are isotopic iff the braids are 
connected by a finite sequence of Markov moves. Given any h £ _B„, a type I Markov move 
consists of conjugating b by any braid x G i?„, taking b t-^ xbx~^. A type II Markov move 
adds another strand to b which is twisted either clockwise or counterclockwise with the n'th 
strand at the bottom of b, acting as b 1-^ a^L{b). It is clear that the trace closures of braids 
related by the first Markov move arc isotopic, since the top braid x is untied by the bottom 
braid x~^ . A picture more clearly illustrates the invariancc of the trace closure under the 
second Markov move, as shown in Figure [5T4l 

These topological considerations are refiected algebraically as follows. It is a theorem of 
Ocneanu [16] (see [19] for a simple proof) that for every 1 <k < I and for each ?i, there is a 
linear function trfc^ : Hn{q) C uniquely determined by 

1. trfc£(l„) = 1 for 1„ G Hn{q) 

2. tTkeiwv) = tvkiivw) for w,v e II„{q) 

3. tTkeieniiw)) = [fc]7[2t trfc^(w) for aU w G for aU n. 
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b b 




Fig. 5.4. Trace closure after the second Markov move on some b G B-j. The dotted box on the 
left represents the braid o"3(,(f)) £ B4. We have omitted the strands which close the two leftmost 
strands for visibility. 



Using (|3.7p . a direct calculation shows that 3. implies that for w G Hn{q), 

trfcf (g„t(u;)) ^ (q - (1 + q) ^^^^2]l ) ^^''^^^^^ = "[fc]7 *i"'^^H- (^-l^) 

Such a function is known as a positive Markov trace on the inductive limit Hao{q) of the 
Hecke algebras Hn{q)- This trace has a representation-theoretic formula which we recall 
below in (|5.2ip . But first, we show how it can be used to express the Jones and one- variable 
HOMFLYPT polynomials. One first extends this trace to S„ via the identification af^ ^ 9^ ■ 
Cyclicity of the trace implies that trfc^ is invariant under type I Markov moves. However, it is 
not invariant under the second move because (|5.17p implies that multiplication of w G Hn{q) 
by either g„ or affects the trace nontrivially. We will therefore construct, for each n, a 
normalized version trfc^ : Hn{q) ^ C of the Markov trace which is also invariant under type 
II Markov moves, i.e. 

Ffe,(a^V(&)) =t?M(6) (5.18) 
for each b G For this, we define 

t^kiib) = [kYr'<l~'^"^'^^Mb)- (5.19) 

Here, e{b) is the sum of the exponents of the generators in the expression of b. Now if 5 G Bn 
and b' = cr,f ^i(^) ^ B^+i-, 

t^kiib') = [ky;q-'i^-^'>'hTki{b') 

= [fc]^,-^(^W±i)^tr..,(6) 
[k\t 

as required. To see that this normalized trace is equal to the one-variable HOMFLYPT of 
the trace closure of 6, we follow [19j by considering an arbitrary oriented link diagram Lq 
described as the trace closure of some braid. Choose a crossing in the braid arbitrarily - 
we assume it to be positive without loss of generality. By a type I Markov move, one can 
then express Lq as the trace closure of a braid where the given positive crossing is at the 
bottom, i.e. Lq = <Tib. Inserting a positive or negative crossing adjacent to the one specified 
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above results in the oriented links L^^ = a^b and L_ = b. Using the formula ()5.19|) to express 
the normalized trace in terms of the linear trace tr^^, and replacing af using the quadratic 
relation p.4p . it is straightforward to check that 

Therefore, the normalized trace satisfies the definition (|2.2p of the one-variable HOMFLYPT 
polynomial from the skein relation (|2.ip with (g) ~ tr^f ((7^6), so that 

Hl'hq) = m^-'q-'^'^'hiMb). (5.20) 
For w G Hn{q), the Markov trace has the following representation-theoretic formula 

tTMiw)^ J2 4''''Tr^f^)(H. (5.21) 

(k £) 

Here, the Markov weights s\ ' are Schur functions given explicitly by 



n lj - t-t f^\g . . 



for each A € A„ ' and h{i,j) denotes the hook length of the box with row-column coordinates 
(i, j) in A. For each A G A^'''^^ if we write d^*"'^' for the dimension of vj^'''^\ we may evaluate 
()5.2ip at the identity 1„ G Hn{q) to obtain the following formula relating the Markov weights 
and the dimensions of the irreps: 



^ 4'='^)df^'=tr,,(l„) = l. (5.23) 



aga: 



We are now equipped to verify the bound on the HOMFLYPT given at the end of Section [2] 
Indeed, by unitarity we may combine (|5.20|1 . (|5.2ip and (|5.23p to obtain, for each b G i?„, 



AGAk''''' 



Jkj),{k,e) 

''X 



AgA 

ITl-l 



= [k]r- (5.24) 

Specializing to fc = 2, we can obtain a more explicit formula for the Jones polynomial. 
The appropriate Markov weights Sj^^ when Ai + A2 = n can be computed as follows. By 
Lemma 3.5(b) of |3H], if ^' is the (2, ^)-diagram obtained by adding r columns of two boxes to 
the left of another (2, ^)-diagram A, it follows that s^^,'^'^ = s'^'^\ A direct computation 

from the hook length formula (j5.22p gives the Markov weight for a single row diagram [m] as 
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Since A = [Ai, A2] can be obtained by adding A2 columns of two boxes to the left of the row 
diagram [Ai — A2], it follows that 



„(2/) 



(2/) _ '^[Ai-Aa] _ [Ai - A2 + l]e _ [Ai — A2 + l]e 



[2] 



2A2 



[2] 



A1+A2 



mi 



Plugging into (|5.19p thus leads to the following expression for the evaluation at g = Q^-m/e q£ 
the Jones polynomial of the trace closure of a braid b E Bn'- 



1 



\e{b) 



J2 [Ai-A2 + l]fTr7rf')(6) 



AGAi; 



(5.25) 



Because the numerator of the coefficient in front is just a complex phase, the absolute value 
of the Jones polynomial is 



|j.(e2"/^)| = |F2,(6)U 



1 



E 

AeA<.='' 



fAi 



A2 + l]fTr7rf')(6) 



(5.26) 



Jones polynomial of generalized closure from the Temperley-Lieb algebra 

Here we obtain representation-theoretic formulae for the Jones polynomial of the generalized 
closure of a braid. For this, we need to briefly describe the connection between the Jones 
polynomial and the Temperley-Lieb algebra. 

The Temperley-Lieb algebra 

As we have described in the previous subsection, the representations corresponding to the 
case k = 2 are all that is relevant for the Jones polynomial. In these cases, it is known that 
the representation matrices satisfy tt^'^'^^ {eiBi+iei — re,) = 0, where r = [2]^^. This can be 
shown using the fact that each of the q-deformed three-row anti-symmetrizers 



1 - <? ^{9i + 9t+i) + q ^(,9i.9*+i + g-i+igi) 



1 9igi+i9i 



vanishes on representations with fewer than three rows (the proof is similar to the well- 
known undeformed case at g = 1, where there can exist no totally antisymmetric state of 
three qubits). The quotient of the Hcckc algebra by the relations e^ei+iei = re^ yields 
the Temperley-Lieb algebra Ti„(T) which is generated by 1 and idcmpotcnts {ei, . . . , e„_i} 
satisfying 



e^ 

c-iej ei 



ei 

TBi 

ejBi 



.?l = l 



(5.27) 
(5.28) 
(5.29) 



Temperley-Lieb algebra is advantageous because it provides representations of the diagram 
monoid Kn [21], defined as follows. For each 1 < i < n — 1, define the "cup-cap" diagram 



LUi 



18 The Jones Polynomial: Quantum Algorithms and Applications in Quantum Complexity Theory 



where the cup-caps act on strands i and i + 1. Multiplication of cup-caps is performed 
pictorially (up to isotopy), amounting to stacking the generators tOi in the same way as with 
the generators of Then, Kn is generated by the cup-cap generators {wi, . . . , ti;„_i}, 

together with an extra closed-loop diagram S, which satisfy the relations 

= Suj, (5.30) 
LOiUjjUii = uii |z-j| = l (5.31) 

ojiUjj = ujjijji |i— j| > 1. (5.32) 

These relations are perhaps best understood in terms of pictures which we invite the reader 
to draw, or otherwise to consult Figure 8 of [H]. While these relations are similar to those 
used above to define the Temperley-Lieb algebra, a more transparent connection is obtained 
by reformulating (|5.27p - (|5.29p in terms of the scaled projections Ei = [2]eei, in which case 
the corresponding relations 

Ef - [2]eE, (5.33) 
EiEjE, ^ E, \i-j\ = l (5.34) 

E,Ej = EjE, \i-3\>l (5.35) 

are formally identical to (|5.30[) - (|5.32[) . The diagram monoid Kn is then represented inside 
TLn^r) via the map taking coi i— > Ei and (5 i— > [2]^. As with _B„, we use this map to obtain 
representations of Kn from those of TLn (t) by setting 

For every even number 2p, define the rectangular tableau t2p G T^p.p] by 



t2n 



1 


3 


5 


2 


4 


6 



2p 



(5.36) 



We slightly abuse notation and write the projection onto the corresponding vector |t2„) as 
t2p = \t2p){t2p\- We pause to mention a convention; given any t G T[pp] t' g T^'^\ let tt' 
be the tableau of shape [p + Xi,p + A2] obtained by adding 2p to each box of t' and placing it 
to the right of t. It is easily checked that tf € t!"^'^} , , This defines a natural embedding 
of vector spaces V[p^p] Va ^ V[p+x^^p+x^] which we denote \t)\t') 1— > \tt). Throughout, we 
freely make this identification on vectors and operators. Now consider the following lemma. 
Lemma 5.1 Let p and n he positive integers satisfying 2p < n. If X = [Ai, A2] G then 

I otherwise. 

Proof: Abbreviating Ci = ttI^'^' (ej, recall that ej\t) ~ \t) iff j and j + 1 arc in the same 
column, while ej\t) = iff j and j + 1 are in the same row. Because the 62^-1 are mutually 
commuting, if 2i — 1 and 2i are in the same row of a given tableau for some \ < i < p^ \t) is 
annihilated by 6163 • • • e2p-i, and thus also by ojiuj^, ■ ■ ■ aj2p-i- Therefore, only vectors \t) for 
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which t ~ t2nt' for some t' G T^x f-p x^-p] survive. Clearly this is impossible if A2 < p. On 
the other hand, for 1 < i < p the 62^-1, and therefore the W2i-ij may only act nontrivially on 
the leftmost p columns of any given tableau. By the convention mentioned above, this means 
that when A2 > p, the quantity in the lemma is proportional to the tensor product of the 
rank 1 projector t2p with an identity matrix as required. □ 



It is possible to combine i?„ and Kn into a unified structure which was called a braid 
monoid by Kauffman |21j . As this braid monoid contains diagrams consisting of both cup- 
caps and twists, its generators are the union of those of i?„ and Kn. However, there are 
considerably more relations that these generators must satisfy (see e.g. [6]) in addition to those 
of the braid group and diagram monoid. For our purposes, however, we will not need the full 
braid monoid. Instead, we shall only need to consider such generalized braids of the form hw^ 
where h e _B„ and w e Kn is a word in the o^.^'s, in which case nn'^^bw) — T:ii'^\b)'Kl?'^\w). 

Jones polynomial of the generalized closures 

For any braid b G -B„, we give a formula for the absolute value of the Jones polynomial of 
the generalized closure of b when a; = y = 1„ e B„. Such a closure then depends only the 
parameters p and r which satisfy 2p + r ~ n. We obtain this formula by computing the 
normalized Markov trace of the braid b with p cup-caps at the tops of the leftmost 2p strands. 
The formula is 



tT2e(bujiUJ3 ■ ■ •a'2p-i)| 

E [Ai - A2 + l]iTT (vrf ■')(&)^f '^(c.ia'g • • -^2^-1) 



1 

M 



[2] 



E [Ai - A2 + 1]^ Tr (tt^ (6) it2p 1 A_ [p,p] )) 



[2] 



[/.I -^l2 + l]i Tr (^g;p]V^(5)(i2p ® 1^)) 



(5.37) 



In the second line, we have used the formula (|5.26p for the absolute value of the normalized 
Markov trace. The third line is by Lemma 1 5. 11 while the last line is straightforward. 

We remark here that this expression yields an immediate upper bound on the absolute 
value of the Jones polynomial of a link obtained as a generalized closure of a braid. By 
unitarity, it is clear that the matrix traces above are maximized when the formula is evaluated 
at the identity braid 1„ G -B„. In this case, we obtain 



Ai2 



MSAi.- 



[2] 



p+r— 1 



(5.38) 
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where we have used (j5.23|) for the last step. By speciahzing (|5.37p to the case where r = 
and p = n, we obtain the foUowing formula for the Jones polynomial of the plat closure of 

b e B2n- 

J{le^^-/')\ ^ [2]ri|Tr(^£3(5)t2„) 

= [2]ri|(i2„k£15(6)|t2„) . (5.39) 

6 Algorithms for approximating the Jones and HOMFLYPT polynomials 

In this section, we use the circuits of Section [4] to give quantum algorithms for approximately 
evaluating, in the sense described in Section^ the Jones and one- variable HOMFLYPT poly- 
nomials at primitive roots of unity. In particular, the results of this section constitute a proof 
of Theorem 12. II We begin by giving a quantum algorithm which obtains an additive approx- 
imation of the one-variable HOMFLYPT polynomials H^'^'f of the trace closure of a braid, 
generalizing the Jones polynomial algorithm of [2]. If fc is a fixed constant, this algorithm 
runs in time polynomial in the length of the braid. Then, we sketch two complementary ways 
in which this algorithm can be generalized to approximate the absolute value of the Jones 
polynomial of any generalized closure of a braid. 

Approximating the HOMFLYPT of the trace closure 

Recall the setup of Problem l2.1l We are given a length-m braid b £ Bn, positive integers fc, £ 
satisfying 2 < k < £, and a number S > 0. We assume that A: is a fixed constant. Because a 
length m braid can act nontrivially on at most 2m strands, we lose no generality in assuming 
that n < 2m. We showed in (jOi)) that the maximal value of | i/i''^ (e^'^^/O | 

over all b G Bn is 

We therefore focus on estimating the following normalized version of the one- variable 
HOMFLYPT which, by ((OO)) and (jOTjl . has the following expression for 6 G S„: 



The Markov weights s\ ' are defined in (|5.22p . We sample from a complex- valued random 
variable with the above expectation by a procedure similar to that of [5] . First observe that by 
(|5.23p . the numbers |Px = s^*"'^^ S^'^^ : A g aI*"'^^! are probabilities. Because k is constant, 

the number of weights grows polynomially with n, so choosing a A G aI^'^^ with probability 
Pa can be done efficiently because those probabilities can be precomputed and sampled from 
by a standard procedure. Second, note that given A € A'k''^\ there is an efficient procedure 
for choosing a path t € x^'^^ uniformly at random. This is because the Young graph A^^^f^ is 

a layered graph, with a polynomial number of diagrams < ( min{n, £ — fc} + l)*" in each layer. 
Therefore, the number of paths from to a given node in the diagram can be computed in 
advance. To choose a path ending in a particular diagram A € A^*^'^', the trick is to move in 
reverse; starting at A, choose the next diagram A' £ A^^^l with a probability proportional to 
the number of paths from to A', divided by the total number of paths from ending in the 
n — I'st layer. 
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(k 

The algorithm proceeds as follows. One first chooses a tableau t e Tn ' randomly as 
above. Next, one samples from random variables X,Y Cz {±1} with conditional expectations 
satisfying 

E[X + iY\X,t] = {t\TT[^-^\b)\t). 

This can be done efficiently using the quantum circuits for the Jones- Wenzl representation 
described in Section |4l together with the following standard lemma which we prove in the 
appendix for convenience: 

Lemma 6.1 (Sampling Lemma) Let U be a quantum circuit of length 0(poly(n)) act- 
ing on n qubits, and let ji/;) be a pure state of the n qubits which can be prepared in time 
0(poly(n)). It is then possible to sample from random variables X,Y £ {il} for which 

E[X + iY] = {ip\U\ilj) 

in 0(poly(n)) time. 

This process is repeated, averaging the values of X and Y , until the desired precision is met. 
The nonconditional expected value oi X + lY is 



A6A 



{k,i) "A 



E[X + iY\X,t] 



ten 



= E 4''^ E (i\-i'-'\b)\t) 



Therefore 



e E[X + iY] 



so a sufficiently good estimate of E[X + iY] provides the needed estimate of the normalized 
one-variable HOMFLYPT polynomial. By the following lemma, we will see that the above 
procedure only needs to be repeated polynomially many times to obtain the desired estimate. 



Lemma 6.2 (Chernoff Bound) Let {Xi, . . . ,Xm} be real-valued random variables satisfy- 
ing \Xj \ < 1 and EXj = ii. Then 



Pr 



1 



> 5) <2 exp 



MS"' 



Suppose the above algorithm is run M times, obtaining random variables {Xj, Yj}jLi- Setting 
Zj = Xj + iYj, the union bound and the above lemma imply 



Pr 



1 



>S} < Pr< 



< 4 exp 



1 ^ 



M6^ 




1 



3 = 1 



> 



V2 
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Fig. 6.1. The trace closure of 6 G B3 is isotopic to the plat closure of the modified braid 
(<t| o-|) ~ ^ fecrl cr| = (T^^o'^'^fT^^6(T3(T4fr2. Note that we identify b with its image t o t o t(6) G Bg. 



Thus, for any > 0, the probabihty that Z = e 21 Zj deviates from the desired 

value of the normahzed one- variable HOMFLYPT polynomial by more than S is smaller than 
the desired error of 1/4, provided that the number of samples satisfies M > 32 ln(2)/(5^. Since 
each sample is obtained in 0(poly(?Ti)) time, the desired additive approximation is obtained 
in 0(poly(m, 1/(5)) time as required. 

Approximating the Jones polynomial of generalized closures 

In this section we show how to approximate the absolute value of the Jones polynomial of 
generalized closures of a braid. As mentioned in Section [21 we focus on the absolute value 
because there is no unambiguous way to assign an orientation to a generalized closure. Recall 
the setup for Problem 12.21 We are given a length-m braid b € two nonncgative integers 
p, r for which 2p + r = n and a positive integer £. We lose no generality in assuming that the 
top and bottom braids which define the generalized closure are x,y ^ !„■ As we showed in 
(|5.38p , the largest absolute value of the Jones polynomial is obtained when the closure yields 
p + r trivial knots. Therefore we focus on estimating the normalized quantity 



p-\-r—l 



[2] 

We show two different ways of achieving this goal, as each constitutes a different generalization 
of the previous algorithm which is interesting in its own right. The first makes use of the 
following proposition [19], which shows that the generalized closure ^^'''(6) can be written as 
the plat closure of a related braid on 2p + 2r strands. 

Proposition 6.3 Let b e i?2p+r- Then ^i'l{b) = c^^bc, where c = o'2p+r'^^^ ' ' ' '^2pt3'^2p+2 
and for i < j , we write aj ~ aiai-^-i ■ ■ ■ <Tj_20'j-i for the braid where strand j is moved above 
its neighboring strands on the left and is inserted in position i. 

While we invite the reader to verify this, the main idea is contained in Figure [^?T1 where 
we illustrate how the trace closure can be obtained from the plat closure. Note that since the 
braid c in the proposition consists of * " *^(^^) generators, the modified braid b' = c^^bc 
is at most polynomially longer than b itself. Setting p' = p + r, observe that since b' G B2p' , 
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(|5.39p implies that 



1 



[2] 



p+r— 1 



In the same manner as with the HOMFLYPT. one then samples from a complex random 
variable Z satisfying EZ = {t2p'\T^^i^p>](b')\t2p') ■ After averaging the results a suitable number 
of times, the absolute value of the average is reported. We omit the details as they are 
identical to the former case. 

Another procedure for this task makes direct use of the following representation-theoretic 
formula for the absolute value of the Jones polynomial: 



1 



J(^f;[(6),e2-/^) 



[2] 



[2] 
1 

w 



p+r-1 



J2 [Ai - A2 + 1], Tr (nl^il+x (b) (hp ® 1 a 



E 



[Ai 



A2 + l],Tr(7rg;5^,(&)(i2p®lA)) 



,(6.1) 



Recall that by (|5.23[) . the numbers |^a = s^^'^^d!\'^^ : A € A^^'^-*! sum to unity and thus define 



probabilities. The algorithm proceeds as in [5] by first selecting an irrep A G A^^'^-* at random 
with probability P\ = s\' d\' . Then, a path t in the Young graph A}^ is chosen uniformly 
at random. By the same arguments as in the HOMFLYPT case, one may then sample from 
random variables X,Y Cz {±1} with conditional expectations satisfying 

E[x + tY\x,t] = {t2pm7Tl^;2+^{b)\t2p)\t). 

This procedure is then repeated, averaging the results until the desired precision is attained, 
after which the absolute value is reported. Using (j6.1|l . it is straightforward to show that 

1 



E[A: + iy] 



Mr-' 



j(^f;[(6),e2W^) 



The remainder of the proof is identical to the HOMFLYPT situation. 
7 Simulating quantum circuits with braids 

In this section, we show how to "compile" any quantum circuit on n qubits consisting of m 
two-qubit gates into a description of a braid on 4n strands with poly(m) crossings in such a 
way that the output of the circuit is encoded into an evaluation of the Jones polynomial of the 
plat closure of the braid at any primitive fth root of unity with £ > 5, £ 6. While this was 
first proved in [12j , our proof is simpler as it only uses the representation of corresponding 
to the rectangular diagram H+H . 

There are uncountably many quantum circuits on n qubits, but countably many braids. 
Therefore, we will have to settle for approximations of quantum circuits. Given unitaries U 
and U' and some < e < 1, we say that U e- approximates U' if \\U — U'\\^ < e, where for 
a square matrix A/. ||Af ||^ denotes the operator norm, or largest singular value of M. Our 
main goal in this section will be to prove the following: 
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Theorem 7.1 Let U ~ UmUm-i ■ ■ - Ui be a quantum circuit of length m acting on n qubits 
and let e = r2(2~P°''''")) . Then there is a braid b G B/^n with 0(poly(m)) crossings such that 
U is e- approximated by '''[|„^2n] (^) • case, the following inequality holds: 



<e. (7.1) 



Moreover, such a braid can be found from a description of the circuit in 0(poly(TO)) time on 
a classical computer. 

To prove Theorem 17.11 we must first introduce some machinery. A finite set of gates in 
U{H) is said to be universal if it generates a dense subset of U{H). The following lemma 
implies that the measurement probabilities of an e-approximation of a quantum circuit differ 
by at most e from those of the desired circuit. The proof is in the appendix. 
Lemma 7.2 Let U and U' be d x d unitaries for which U' e- approximates U , where e < 1. 
Then, for every pair of normalized pure states \4>) G C'', we have 

\Mum''-mu'H)?\<e- 

The next theorem assures that universal sets of unitaries suffice to efficiently e-approximate 
an arbitrary quantum circuit. 

Theorem 7.3 (Solovay-Kitaev) Let Q C U{d) be a finite universal set of unitaries on a 
Hilbert space of constant dimension d which are closed under inverses, i.e. Q^^ = Q. Then 
any unitary U G U{d) can be e- approximated by a sequence o/ 0(poly log(l/e)) gates in Q. 
Moreover, such a sequence can be found by a classical computer in 0(poly log(l/e)) time. 

Fixing any finite, universal set of two-qubit unitaries G2 C U{4), let Qn C C/(2") be the 
set of gates obtained by letting Q2 act on all pairs of neighboring qubits. Given a length m 
circuit UmUm-i ■ ■ - Ui consisting of gates from C/„ and any fixed desired degree of accuracy e, 
the Solovay-Kitaev Theorem further implies that for each i, there is an e/m approximation 
U[ to Ui which is obtained by a sequence of 0(poly log(m)) gates in Q, and also that such 
a sequence can be found in 0(poly log(7T!,)) time (since e does not grow with to, it does not 
figure into the complexity estimate). Observe that for any unitaries U and U' , the operator 
norm is stable, meaning that ||C^ — C^'foo = \U ®^d — U'®\d\^ for every finite dimension d. 
Therefore, the same approximations hold for the embeddings of the local gates into C/(2"). 
By a telescoping sum and the triangle inequality, we have 

||C/™ • • • C/2C/1 - ■ ■ ■ U!,U[ lU < E i^' - ^^ioo- (7-2) 

1=1 

Therefore, there is an e-approximation to the desired circuit which uses 0(TOpolylog(77i)) 
gates from Q^. So we see that working with finite universal sets of gates is just as good as 
working with all possible two-qubit unitaries, provided that one can live with bounded errors 
and polylogarithmically more gates. In Section [8l we define classes of quantum computational 
problems which are insensitive to these limitations. 

Throughout this section, we fix an £ > 5,^ 7^ 6 and for each diagram A, we abbreviate 
tta = TTj^ ' ' and Vx = . We begin by describing how to encode a single logical qubit into 
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Fig. 7.1. Qubit path basis for Wi. The upper path corresponds to the basis state jio), the lower 
path to 



Vg. This is a two-dimensional representation of B4 with orthonormal basis 



1 


2 


3 


4 



These basis states may also be viewed as paths in the Young graph (see Figure 17. ip . Let t 
and t' be (2,.£)-tableaux of respective shapes and [/c, k]. Recall the convention mentioned 
before Lemma ISTT] where we write tt' for the tableau of shape [j + k,j + k] obtained by adding 
2j to each box in t' then placing it to the right of t. Central to our proof is that tt' will 
always be a (2, ^)-tableaux. Recall that this allows us to define a natural embedding of vector 
spaces Vjj j] (8) Vj^ ^j ^ Vy+kj^/.] via the mapping \t)\t') 1— > \tt'). We encode n qubits into a 
subspace T/F*^"-' of the irreducible representation V\2n,2n\ of which is spanned by the 2" 
computational basis states Wt^^tx^i ■■■'tx„) ■ {xi,X2, ■ ■ ■ ,Xn) S {0,1}"}. By our embedding 
the subspace W^'^") can be identified with a tensor product (^"^j^ Wi, where for each 1 < i < n, 
we identify Wi with V^, so that it is spanned by \to) and In this subspace of V[2„,2n]i 
each consecutive group of four strands is associated to a single qubit in the following way. 



For each l<j<A:<n— 1 



define the subset Y^j 



, cTfe} of generators of B„ 



and write C i?„ for the subgroup generated by S*^ ^. For each 1 < i < 71 — 2, the image 



'J — ^•^'^ yj^^^ ^^.^i-^^^^^ ^^^^^^^^^^ 

of the 8 strand braid group under the representation Tr[2n.2n]: when restricted to 

acts nontrivially only on the subsystem Wi (g) Wi+i. In fact, it rapidly approximates 
any U € SU{Wi (8) Wi+i). To see this, it suffices to look at the case i = 1. Note that 
Wi ® W2 C Vrm- The embedding (g) ^ %m looks like 



IMIto) ^ 

\to)\ti) ^ 

\tl)\to) ^ 

ltl)ltl) ^ 
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While this proves that dim Vfm > 4, we mention that its dimension is 13 when £ = 5 and is 
14 when £ > 6. So it is clear that SU{Wi (g) W2) C SUjVfpm). Consider now the following 
proposition. 

Proposition 7.4 (Freedman, Larsen, Wang [13]) Let £ > 5,£ ^ 6. The image of Bg 
under tt^^ is dense in SU (v^^^'^ ^ . 

With this density result in hand, we are now equipped to prove Theorem 1 7. II 

Proof of Theorem I7.lt Set e' = e/m = fi(2^™/m). Together with the Solovay-Kitaev 
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Theorem, Proposition 17.41 implies that for each U' G SU{Wi (E) W2), there is a braid b G Bs 
of length 0(poly log(l/e')) = 0(poly(m)) for which nfm(b) e'-approximates U' , where we 
consider U' to be embedded in SU (Vprn). Note that while some unitaries Tr^{ai) will not 
be contained in SU{Wi W2), the construction above ensures that Trnm(b) can be made 
arbitrarily close to any U' G SU{Wi (g) M^2)- 

Because of the four-periodicity of our encoding, the same arguments apply for approxi- 
mating gates on any other adjacent pair of qubits. Now, let U = UmUm-i ■ ■ - Ui as in the 
statement of the theorem. Without loss of generality, we may assume that each Ui has unit 
determinant. By (|7.2p . we may conclude that there is a braid b E i?4„ of length 0(poly(m)) 
for which TT[2n.2n]ib) e-approximates the circuit U, and also that such a braid can be deter- 
mined on a classical computer in 0(poly(m)) time. The rest of the theorem now follows by 
Lemma 17.21 and the expression (|5.39p for the absolute value of the Jones polynomial of the 
plat closure; the normalization constant l/[2]^"~^ is so because the braid has An strands. □ 



8 Complexity-theoretic applications of the Jones polynomial 

Complexity theory characterizes the asymptotic consumption of resources, such as time and 
space, required to obtain the solutions of certain classes of problems on a computer. Writing 
{0, 1}* for the set of all binary strings of finite length, a decision problem asks for the evaluation 
of some function /: {0,1}* {0,1} on an arbitrary input. Such a function determines a 
language L C {0,1}* defined as L = /^^(l), so equivalently, the task is to determine, for 
each bit string x", whether G L. Examples of languages include the set of all binary 
strings of even parity and the set of prime numbers, written in binary. A complexity class is 
a collection of computational problems, usually defined in terms of the resources required to 
determine membership of arbitrarily long bit strings in each language contained in that class. 
For instance, any class of languages is a complexity class. Two well-known complexity classes 
are P and NP. The former is defined to contain every language L C {0,1}* for which there is 
a uniform family {C„ : {0, 1}" {0, 1}} of polynomial-size classical circuits (see e.g. [501 15T] 
for a definition of a classical circuit) for which C„(a;") = 1 when cc" € i, while C„(x") = 
otherwise. By a uniform family, it is meant that there should be a classical algorithm which, 
given the number n as input, produces a description of the circuit C„ in 0(poly(n))-time. 
We mention that this definition of a uniform family of circuits extends to quantum circuits 
in the obvious way. NP, on the other hand, contains those languages L for which there is a 
uniform family of classical circuits {C„ : {0, 1}" x {0, 1}™ {0, 1}}, where the number m of 
guess bits grows polynomially with n, so that for each x" e L, there exists a witness string, 
or certificate y™, so that C„(x", y™) = 1 if x" € L, while if x" ^ L, Cn{x",y"^) = for every 
possible certificate y"\ Below, we will recall the definitions of certain quantum generalizations 
of these classes, showing how they may be characterized using the Jones polynomial. 

The complexity class BQP is defined to contain those decision problems whose solution 
can be determined with bounded error in polynomial time in the standard quantum circuit 
model. Formally, a language L is in BQP if there exists a uniform family of quantum circuits 
{[/„} of polynomial length such that 



(a;"00 



■•0|[/„|x"00...0)p|^J/^ 



if x" e i 

if x" ^ L. 
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This relates the usual definition of BQP, in which a single output qubit is measured to 
obtain the result of a computation, as follows. Consider a circuit such that if x" G L, 
a measurement of the first qubit of U^\x"0---0) in the computational basis yields 1 with 
probability > 3/4 while if ^ L, that measurement outcome will be with probability 
> 3/4. Appending an additional qubit 10)"* and labeling the output qubit of as B, we 
then take 



I' 



\i)a\ 



'){i^®u':). 



Note that the Solovay-Kitaev theorem implies that the class BQP remains unchanged if 
the particular universal gate set chosen for the model merely contains a universal set of gates 
for each pair of qubits. The following theorem essentially appears in [7j; we prove it in this 
paper for the sake of completeness. 

Theorem 8.0 Let A be a random oracle which solves Approximate Jones Closure with 
(5 = 1/8 for the plat closure and some con stant £ > 5,^ 7^ 6. Then P^^ = BQP. Loosely 
speaking, this means that Approximate Jones Closure is "BQP -complete. " 
Proof: To begin, it is an immediate consequence of Theorem 12.11 that P"* C BQP. To see 
that BQP C P'^, let L G BQP and let U be n'th quantum circuit in the uniform family of 
circuits for L. By definition, U is of length m = 0(poly(n)) and acts on n' = 0(poly(n)) 
qubits. For each x", define the circuit Wx^ = ^^=1^^', which acts as the identity on the 
ancillary qubits. Because 

{00■■■0\Wx'^UWx^\00■■■0) = (a;"00---0|;7|x"00---0), 

it follows by Theorem l7.1l that for any constant e > 0, there is a braid b of length 0(poly(m)) = 
0(poly(?i)) on 4n' strands which satisfies 



J{b, e 



2Txill\ 



|(a;"00---0|C/|a;"00---0)p 



< 



Moreover, there is an 0(poly(n))-time classical algorithm for determining such a braid from a 
description of W^^ UW^^ ■ By invoking the oracle A, we obtain a random variable Z satisfying 



Pr 



[2] 



2n' 



J(b, e 



< 1/ 



> 3/4. 



Since — b^\ < 2\a — b\ for < a, 6 < 1, the last two estimates can be combined via the 
triangle inequality to show that | - | (a;"00 • • • 0|[/|a;"00 • • • 0) p | < 1/4 + e with probability 
at least 3/4. Choosing e = 1/8 is sufficient to complete the proof. □ 



The complexity class QCMA is a certain quantum analog of NP, consisting of those lan- 
guages L for which a computationally unbounded oracle (generally personified as Merlin) can, 
with high probability, efficiently convince a computationally bounded verifier (Arthur) that 
each G L is actually contained in L. Cheating is also discouraged, in that if x" ^ L, 
the probability that Merlin will convince Arthur otherwise should be small. Formally, we 
say that a language L is in QCMA if there is a uniform family of quantum circuits {Un} 
of polynomial length such that for each length n string .t" e L, there is a length m binary 
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String for which P(a;",y"') = |(a;"y™00 • • •0|C/„|x"?/™00 • • -0)12 > 3/4, while if ^ L, 
we have y™) < 1/4 for every such y™. Wc remark that the usual definition of QCMA 

(see e.g. [3] for the original definition) is phrased in terms of measurement probabilities of a 
single output qubit; by the remarks following the definition above of BQP, the definition we 
give here is equivalent. For the same reasons as with BQP, the exact value 3/4 of the success 
probability is not important; what is important is that the square of the matrix element is 
strictly bounded away from 1/2. 

Given a complexity class C, a problem is said to be C-hard if it is at least as hard as 
any other problem in the class C, in the sense that access to an oracle which immediately 
computes the solution of that problem allows the solution of any problem in the class C in 
polynomial-time. If that problem is also contained in the class C, we say that it is C-complete. 
We will show that the following problem is QCMA-complete. 

Problem 8.1 (Increase Jones Plat) Given a braid b E B2n with 0(poly(?i)) crossings, 
a fixed integer i > 5, £ ^ 6, and a class of braids C2n C i?2ra of length 0(poly(n)) for 
which membership can be decided in 0(poly(n)) time on a classical computer, decide, with 
the promise that only these two cases can occur, whether there exists another braid c € C2n 
for which the absolute value of the normalized Jones polynomial of the plat closure 

^ J(^*?i,e2"*/0 ' > 3/4 



[2]," 

or if it is < 1/4 for all braids c G C2n ■ 

Theorem 8.1 Increase Jones Plat is QCMA-complete. 

Proof: To begin, wc demonstrate that Increase Jones Plat is in QCMA. Given a braid 
& € B2n of length 0(poly(ri)) , Merlin begins by sending Arthur the classical description of the 
"witness" braid c. Next. Arthur checks to sec if c G C2n, which can be done in polynomial-time 
by definition. Because of the promise, Arthur only needs to learn the value of 

2 



[2]"- 



with an accuracy of < 1/2 with probability < 3/4. For this, it is sufficient for him to 
run Approximate Jones Closure on a quantum computer with 5 = 1/4, which takes 
0(poly(n))-time. 

To see that this problem is QCMA-hard, and is thus QCMA-complete, suppose that L G 
QCMA. Let U be the n'th circuit from the uniform family of quantum verifiers for L, and 
suppose that U acts on n' qubits while accepting witnesses of size m. Letting e > be a 
constant to be determined at the end of the proof, Arthur precomputes a set of four braids 
Cs C -Bs which e/m-approximate the circuits ^X^^ X^^ : yi,y2 G {0, 1}} (obviously, there 
is nothing to compute when yi ~ y2 = 0). This can be done in polynomial-time. Define 
Cin' C Bin' to consist of braids which act as the identity everywhere except for the strands 
4n 4- 1, . . . 4n -t- 4m, where only a single braid from Cg can act on each group of eight strands 
8i -I- 1, . . . , 8i -I- 8 within those strands. If we define Vy7n in an analogous manner to Wx" in 
the proof of Theorcm l8.01 where it acts only by flipping the witness qubits. it follows by (|7.2p 
that for each braid c in C^n' ■, the corresponding unitary '^^2n'^2n'] ('^) e-a-pproximates Vym for 
one and only one y'" (this is because \\Vym — Vz'^j^ = 2 whenever y™ ^ z'"). Artlmr then 
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tells to Merlin, who sends him a witness braid c from C^n'- Define W^^ as in the proof 
of Theorem 18.01 except that it acts as the identity on the witness and ancillary qubits. By 
the first half of Theorem 17. 11 Arthur may compute a braid b for which the circuit Wx^UWx-^ 
is e-approximated by 7i'[|„^''2n'] (^)- Therefore, we obtain by (|7.2p that ''^[2n'^2n']('^~^^'^) ^ 
3e-approximation of Wj,™ UWx'^ Vy^ . Since 

{m---Q\Wxr^Vyr.UWxr^Vyr.\m ■■■{)) = (.t" ™ 00 • • • | | y " 00 • • • O) , 

we may apply the rest of Theorem 17.11 to obtain that 



(x"y"00---0|J7|a;"2/'"00---0) - 



[2]2«' 



J(c ifoc, e 



< 3e. 



In order to complete the proof, notice that because of the promise, it is sufficient to choose 
e = 1/8. □ 

A rather large complexity class is PSPACE, which contains all languages which can be decided 
using only polynomial space. In [IT], PSPACE is characterized as the class of languages 
which can be decided by applying the same polynomial-size circuit (possibly exponentially) 
many times. Formally, it is shown there that if L € PSPACE, then there exists a uniform 
family of polynomial-size quantum circuits {[/„}, together with a sequence of polynomial-time 
computable natural numbers e„, for which, if /: {0,1}* {0,1} is such that /^^(l) = L, 
then 

(x" 100 • • • 0|C/^" |x"100 • • • 0) = fix''). 

Using this characterization, we will show that the following problem is PSPACE-complete. 
Problem 8.2 (Approximate Concatenated Jones Plat) Given a braid b g i?2n with 
0(poly(n)) crossings, a constant I satisfying ( > 5, £ ^ 6, and a positive integer e which 
is describable by 0(poly(n)) bits, decide, with the promise that only these two cases can oc- 
cur, whether the squared absolute value of the normalized Jones polynomial of the plat closure 
of the concatenated braid b'^ satisfies 



1 



[2]r 



J(6%e'"/') >3/4 



or if it is < 1/4. 

Theorem 8.2 Approximate Concatenated Jones Plat is PSPACE-comp/eie. 
Proof: To sec that this problem is PSPACE-hard, choose any language L G PSPACE, let 
U be the n'th quantum circuit from its associated imiform family and let e = (9(2P°^y(")) be 
the corresponding exponent. Suppose that U acts on n' qubits. For a given x", we will show 
how to use an oracle for Approximate Concatenated Jones Plat to determine whether 
or not x" G L. Define the circuit W^^ as in the proof of Theorem l8.0l except that it also flips 
the n + 1st qubit. By choosing e = e'/e = 17(2"?°'^'^"^) for a constant e' to be determined 
later, we may obtain, via Theorem 17.11 a braid b with only 0(poly(n)) crossings such that 
^[2n^^2n'] E-approximatcs Wx^UW^r^. Note that 



/(x") = (a-"100---0|(W:,,.C/W^«)'=|a;"100---0) = (00 • • • 0|W:,n |00 • • • 0). 
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Because e'/e = e, it follows by (|7.2p that T^[2n'' 2n']i^'^) e'-approximates the cireuit Wx^W^Wx"- ■ 
Choosing any e' < 1/2 is sufiieient to complete the proof of PSPACE- hardness. All that is 
left to do is to check that this problem is contained in PSPACE as well. First, it is clear that 
Approximate Concatenated Jones Plat is contained in BQPSPACE, as our local qubit 
implementation of Section 2] requires only polynomial space. On the other hand, Watrous has 
shown [37] that PSPACE = BQPSPACE, completing the proof. □ 

Another class of computational tasks is comprised of counting problems, which ask, for a 
given function /: {0, 1}* {0, 1}, "how many x" € {0, 1}" satisfy /(x") = 1?" The class #P 
contains those counting problems arising from a function / for which /^^(l) e P. We refer 
the reader seeking more detail to [30j for a more rigorous definition. 

Theorem 8.3 Given a braid b S B2n with 0(poly(n)) crossings and a fixed integer i > 5, 
£ ^ 6, computing the n most significant bits of the absolute value of the normalized Jones 
polynomial of the plat closure of b 



J(6, e^"/^) 



is #T'-hard. 

Proof: Suppose /: {0, 1}" {0, 1} is computable in polynomial-time on a classical com- 
puter. It is well-known (see e.g. [31]), that there is a quantum circuit U, acting on n' = 
0(poly(n)) qubits, which computes / reversibly and exactly, namely 

f/|a;")|y)|00---0) = © /(x") © 1)|00 • • • 0). 

The state 



1= J2 k">|i)|oo---o) 



/On 

a;"G{0,l}'* 



can clearly be prepared efficiently. Observe that 



El 

2" ' 



where Ni = |/^^(1)| is the solution to the #P problem of counting the solutions to /. Let 
W be the unitary circuit which acts with Hadamard gates on the first n qubits, flips the 
next qubit, while leaving the rest alone, and set e = 4^"^^. Clearly V1^|00 • • • 0) = \tp). Since 
|a — 6| < -^/ja^ — 6^1 for < a,b < 1, it follows by Theorem 1 7 . 1 1 that . given a description of 
the circuit W^^UW, we may, in polynomial-time, compute a braid b for which 



I (00 • • • o\w-'uw\oo • • • 0) I - Jib, e^w^; 



In particular, this means that the n most significant bits of the two terms in the above absolute 
value agree. Therefore, an oracle which gives the An' most significant bits of that evaluation 
of the Jones polynomial can be used to obtain all n bits of the number iVi . □ 



In fact, n bits is more than enough bits for #P-hardness. Learning just ^l{n^) of the most 
significant bits, for any constant (5 > 0, is T^P-hard as well. This can be seen by padding with 
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a sufficiently large, but polynomial, number of extra qubits. On the other hand, we will see 
that learning just the highest order bit of the Jones polynomial is PP-hard. 

The complexity class PP, introduced in [T3] , contains those languages recognized by some 
uniform family of polynomial-size probabilistic classical circuits with probability greater than 
1/2. Equivalently, L is in PP if there is a nondeterministic polynomial-time machine which 
accepts G L on more than half of its paths. 

Theorem 8.4 Given a braid b £ B2n with 0(poly(7i)) crossings and a fixed positive integer 
£ > 5, ^ 7^ 6, computing the most significant bit of the absolute value of the normalized Jones 
polynomial of the plat closure of b 



is PP-hard. 

Proof: The proof follows by exactly the same manner as in the proof of #P-hardness. The 
only difference is that the most significant bit of the Jones polynomial is sufficient for deciding 
if the number is solutions is greater than 2""^. □ 

9 Conclusions 

We showed how to implement the unitary Jones- Wenzl representations of the braid group at a 
primitive £'th root of unity in polynomial-time on a quantum computer. We constructed the 
image of each generator using a number of local two-qubit unitaries which was polynomial in 
the number n of strands, independent of how £ grows with n. We then used this model to give 
algorithms for obtaining additive approximations of one-variable HOMFLYPT polynomials 
i/^'^-' on a quantum computer. These algorithms run in time which is polynomial in the 
number of crossings and the inverse of the desired accuracy, independent of how fast £ grows 
with these parameters, provided that k is a constant. While our circuits implementing the 
Jones- Wenzl representations tt'^'^'^^ are of length 0(poly(n, /c)) , it appears that the classical 
part of our algorithm causes the running time to scale exponentially with k. 

Other authors have approached the problem of implementing the Jones representations 
on a quantum computer. In [22], a representation of the three-strand braid group was given 
in a qubit architecture. Another approach toward implementing the Jones representations 
is contained in |35j : however, the final conclusions of that paper rested on an incorrect as- 
sumption that unitary matrix elements could be measured precisely in a single time step. 
Finally, a recent article [14j drew on a connection between topological quantum field theories 
and spin networks in order to outline a quantum algorithm for approximating the colored 
Jones polynomial; those authors left as an open question whether their algorithm could be 
efficiently implemented on a local quantum computer. 

On a classical computer, if the number of strands is constant (i.e. not included in the 
complexity estimate), then the Jones polynomial of the trace closure of a braid with m 
crossings can be computed exactly by simply multiplying representation matrices of the braid 
group, requiring time which is polynomial in m. However, accounting for the number of 
strands, the matrices will be exponentially large in the number n of strands. On a quantum 
computer, however, multiplication of exponentially large matrices is not a problem, provided 
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that they are sufficiently "sparse," or rather, that they arise via local unitary transformations 
of the state space of the quantum computer. On the other hand, there is a sense in which such 
large unitary matrices are not explicitly observable; an exponential number of queries seem 
to be required in order to learn the trace of the matrix representing the unitary evolution 
of a quantum circuit. Nonetheless, the additive approximation of the Jones polynomial we 
achieved with a quantum algorithm is unobtainable on a classical computer, unless of course 
BPP ^ BQP. 

On the other side of the fence, we have simplified the proof from [T^l [H] which shows 
how to simulate the standard quantum circuit model with braids. By using four strands per 
qubit, rather than the three of the original proof, we simplified the representation theory 
dramatically. The four-periodicity of our encoding means that shifting the eight-strand braid 
which represents a given local unitary on two qubits by a multiple of four strands realizes 
the same gate on another pair of qubits. The original result of [12] requires consideration of 
two different cases, depending on whether the first qubit is located at an even or an odd site. 
Our encoding applies to schemes for topological quantum computation based on the braiding 
of nonabelian anyons with label ^ in SU{2)i-2 Chern-Simons field theories with ^ > 5, 
^ ^ 6. On the other hand, our local qubit implementation of the Jones- Wenzl representations 
allows efficient simulation of the braiding of particles with label i in SU {k)i-k Chern-Simons 
theories. While slightly less general than the results of [10], as the only mapping class group 
our methods apply to is the braid group, we have tried to make our presentation accessible 
to readers who are unfamiliar with the intricacies of topological quantum field theories. 
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Appendix A 

Proof of Lemma I6.lt Introduce an extra system labeled C which holds a single control 
qubit and denote the Hilbert space of the qubits acted on by U by A. Because U consists of 
0(poly(n)) gates, then so does the controlled unitary V: AC AC, defined by 1/|-0)|O) = 

|?/')|0) and IV')]!) = Initiahze the qubit in the state "^(|0) + and prepare 

the state Applying the controlled unitary places everything into the state 

|*)-^'^ = -i=(|V')|0) + ([/|v))|i)). 

Writing '0 = | ■(/;)('(/; | , the reduced density matrix t'^ of the extra qubit is then equal to 

= 1 f ^ \ 

^ If 1 

= i(l2 + ajfRe(^|f/|V') +CTyIm(^|C/|V')). 

Since the ±1 eigenstates of ax arc |±) = ^(lO) ± measuring ax on the qubit yields a 
classical random variable AT e { + 1, —1} for whichEX = 'Re{ip\U\4>) as required. Similarly, by 
running the algorithm a second time, we can sample from Y, which satisfies EY = lm{tjj\U\t/j). 
□ 

Proof of Lemma 17. 2t By the triangle inequality, 

\{m^)\' = |(0|t/'|V'> + (0l(t^-c^')IV^>r< |(0|c^'IV'>|' + |(0l(c/-c/')IV'>|'. 

The operator norm satisfies = max{|TrMX| : \\X\\^ = 1}, where the trace norm \\X\\-^ 

of a square matrix X is defined to be the sum of its singular values. Thus, 

\{m~U'm\^ = \Tr{U-U'm{cl>\\' <\\U~U'\\l<e'-<e, 

I 1 2 I 1 2 

Therefore, (t/)|J7|'!/;) — (^IJ/'I-!/;) < e. A similar argument gives the same bound with U 
and U' swapped, proving the lemma. □ 



